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We compute the entropy of a closed bounded region of space for pure 3d Riemannian gravity 
formulated as a topological BF theory for the gauge group SU(2) and show its holographic behavior. 
More precisely, we consider a fixed graph embedded in space and study the flat connection spin 
network state without and with particle-like topological defects. We regularize and compute exactly 
the entanglement for a bipartite splitting of the graph and show it scales at leading order with 
the number of vertices on the boundary (or equivalently with the number of loops crossing the 
^"Y boundary). More generally these results apply to BF theory with any compact gauge group in any 

space-time dimension. 
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Introduction 



o . 

Entropy is a key notion in the search for a quantum gravity, both as the thermodynamic quantity useful to probe the 
physics and potential phenomenology of the theory and as the measure of information useful to identify the physical 
bi), degrees of freedom and their dynamics. Research has focused on the particular case of black holes and has lead to 
the holographic principle directly relating geometric quantities (the area) to the entropy in quantum gravity. 

In the context of Loop Quantum Gravity (see 0] for reviews), most of the black hole entropy calculations have 

■ been performed in the framework of isolated horizons following the seminal work by Ashtekar, Baez and Krasnov 
\Q ' Q. Assuming precise boundary conditions for the gravitational fields on the horizon, they count the number of 
C*") • (kinematical) boundary states consistent with fixing the value of the area. Here we would like to get rid of the 

classical boundary: instead of considering a spin network state with a boundary specified classically, we look at 
a closed spin network and an arbitrary bipartite splitting into inside/outside regions (see e.g |3|). The aim is to 

■ compute the entanglement between these two parts of the spin network for a physical quantum geometry state solving 
<0 the Hamiltonian constraint. The first step that we take here is to work this out for BF theory instead of gravity. Indeed 

BF theories are topological field theories lacking local degrees of freedom, thus allowing for an exact quantization (sec 
e.g [l|). In this context, we know the physical states solving all the constraints and can compute the entanglement 
explicitly. This turns out to be similar to ground state entanglement calculations in some spin models developed for 
topological quantum computation 0, 0] . 

The motivation to analyze BF theory is that it is very close to gravity. First, in three space-time dimensions, 
3d gravity is actually a topological BF theory with the Lorentz group as gauge group. Second, in four space-time 
dimensions, gravity can be reformulated as an constrained BF theory and we can work on a quantization scheme with 
quantum BF theory as the starting point. Of course, the fact that gravity has local degrees of freedom should matter 
in the end. Nevertheless, studying BF theory should allow us to develop mathematical tools and procedures later 
useful for loop gravity. 

In the present paper, we start with a quick overview of BF theory and the definition of the physical quantum states 
as spin network states for the flat connection with possibly particle-like topological defects. Then focusing on SU(2) 
BF theory, we explicitly compute the entanglement between the two parts of such a flat spin network states and we 
show its holographic behavior: it scales with the size of the boundary (more precisely, with the number of boundary 
vertices). We also study the influence of topological defects. We show that they do only affect the entropy when 
located on the boundary between the two regions and we compute the finite variation of entanglement that they 
create. 
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I. AN OVERVIEW OF BF THEORY 

A. SU(2) BF theory: spin networks and physical states 

BF theory is a class of topological gauge field theories defined on a oriented smooth n-dimcnsional manifold M. by 
the following action (see e.g. [l[ for a review): 

S[A,B] = f tv(BAF[A]). (1) 

The gauge group is a (semi-simple compact) Lie group G whose Lie algebra q is equipped with an invariant (non- 
dcgcncratc) bilinear form tr(-, ■). Picking a local trivialization of the principal G-bundle over A4, the basic fields are a 
q- valued connection 1-form A with curvature 2-form F[A] and a g-valued (n-2)-form B. The action is invariant under 
the action of the gauge group, for h E G: 

B^hBh- 1 , A -> hAh- 1 + hdhT 1 . (2) 

It is also invariant under shifts of the B field by an arbitrary (n-3)-form <f>: 

B^ B + d A ^>, 5 A = 0. (3) 

This symmetry kills all local degree of freedom making BF theory a topological field theory. Its classical field equation 
impose a flat connection. F[A] = 0, and a vanishing "torsion" . d A B = 0. In particular, in n = 3 space-time dimensions, 
BF theory for the gauge groups G = SU(2) and G = SU(1, 1) is equivalent to 3d Riemannian and Lorentzian gravity 
(in its first order formulation). 

For our present purpose, we are interested into BF theory from the point of view of its canonical quantization a 
la loop quantum gravity. In the following, we will work with n = 3 and G = SU(2) although all the formalism and 
techniques apply to arbitrary space-time dimensions and arbitrary (compact semi-simple) Lie groups. We perform a 
2+1 splitting of the 3d space-time viewing the 3d manifold M. ~ E x R as a two-dimensional space E evolving in time. 
Then the spatial parts of the connection A and the field B on E are canonically conjugate variables. The theory is 
completely constrained and the Hamiltonian vanishes on-shell. There are two sets of constraints, imposed by the time 
components of A and B, respectively, which turn out to be simply Lagrange multipliers. The first constraints impose 
the flatness of the connection on E. The second set of constraints imposes the vanishing of the torsion on E and is 
usually called the Gauss law. These are first class constraints, respectively generating the translational symmetry ([3]) 
and the gauge invariance ([2]) under the action of the group G. For more details on the canonical analysis and resulting 
structures, the interested reader can check Q . 

The loop quantization scheme is based on a specific choice of wave functions. We choose cylindrical functionals of 
the connection A on E. More precisely, they depend on A through only a finite number of variables: they are functions 
of the holonomies of A along the edges of some arbitrary (finite) oriented graph in E. Considering a particular graph 
r with E and V vertices, we consider the holonomies <7i[A|, ..,<?b[A| of the connection along the edges e = 1, ..,E of 
r and build wave functions of the following type: 

MA) = ^( 9l [A],..,g E [A}). 

Further we require these functionals to be gauge-invariant. Since the action of the group on the connection translates 
into a group action at the end points of the holonomies, g e [A] — > hJ^g e [A]h t ^, where s(e) and t(e) are respectively 
the source and target vertices of the edge e, the gauge invariance of the wave functions tpr involves a invariance under 
the group action at every vertex of the graph T: 

Vh v eG xV , ip(gi, ..,g E ) = ^(K(i)9i h t(i), ■■, K(E)9Eh t (E))- (4) 

Then we need to impose the flatness condition F = on these wave functions. To keep the discussion as simple as 
possible, we consider a graph T which forms a lattice faithfully representing the canonical surface E. More precisely, 
if we cut the surface E along the embedded graph T, then we are left with surfaces all homomorphic to the unit 
disk. These are the faces of the lattice. We now impose that the holonomy of the connection around each face is the 
identity. Thus the projection onto physical states satisfying the flatness condition is implemented by the multiplication 
by 5-function around each face: 



C \eeC ) 



(5) 
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where C labels the loops around the faces of the lattice. This imposes trivial holonomies around every contractiblc 
loop in E while allowing for arbitrary holonomies around the non-contractible cycles of the surface. 

When S is the two-sphere, this gives a single flat quantum state. When S is an orientable compact surface of 
genus n, the physical space of flat quantum states is isomorphic to the space of L 2 functions of 2n group elements 
A\, Bi, .., A n: B n satisfying AiBiA^ 1 B^ 1 ..A n B n A~ 1 B~ 1 = 1 and invariant under diagonal conjugation Ai,Bi i— > 
hAih~ x , hBihr 1 for h G G. For instance, for the two-torus, physical states will be gauge invariant functions of two 
group elements A and B satisfying ABA~ 1 B~ 1 = 1. 

Finally, in 2+1 spacetime dimensions particles create conical singularities, which leave the spacetime flat except 
along their worldline. They are represented as topological defects, which translate into non-trivial holonomies. For a 
spinlcss particle on a given face of the lattice, we replace the (5-function by a ^-function imposing that the holonomy 
around the face has a class angle 9. This deficit angle 9 then defines the mass of the particle (e.g. see 0,13 for more 
details). 

Our goal in the present work is to compute the entanglement on a physical state between a region of £ and the 
rest of the surface and to check whether it satisfies to an "area-entropy" law. On a given graph T, we will consider a 
connected region of the graph and compute the von Neumann entropy of the reduced density matrix obtained after 
tracing out all the holonomies outside that region from a given physical state. In the next part, we will show how 
this works for a BF theory based on a discrete group, in which case the theory can be reformulated as a spin system. 

To conclude this review section, we give a few mathematical details on the SU(2) group. We parameterize group 
elements g as 

g = e ims = co86t + isin9u.a, (6) 

where 9 <G [0, 27r] is the class angle (or half of the rotation angle), u is a unit vector on the 2-sphere indicating the 
axis of the rotation and a are the standard Pauli matrices. Let us point out g(9, u) = g{—9, —ii). In these variables 
the normalized Haar measure dg on SU(2) reads as 



/ dgf(g)= ^ P sin 2 9 d9 [ d' 

JSV(2) Z7T JO JS 2 



uf(e,u). (7) 



By the Peter- Weyl theorem, every function which is invariant under conjugation can be decomposed on the characters 
of the irreducible (spin) representations of SU(2). Such representations are labeled by a half-integer j G N/2 and the 
corresponding characters are 

XM = tr, {pi is)) = Sin( ^+ 1)9 = U 2j (cos6), (8) 
where U n is the n-th Tchebychev polynomial of the second kind. The (^-distribution decomposes as: 

6(g) = ]T (2j + l) Xj (g). 

jeN/2 

Finally, we introduce the distributions Sg that localize group elements on a specific equivalence class under conjugation 
by fixing their rotation angle: 

dgS e (g)f(g) = ^- [ d 2 uf((9,u)). (9) 



-in 
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Its decomposition into characters reads: 



b(g)= E xMxM- (io) 

j'eN/2 



B. A word on Kitaev's spin system and Z2 BF theory 

In this subsection, we give a quick overview of Kitaev's spin system introduced in Q and the corresponding ground 
state entanglement calculations presented in Q. As in the previous section, we consider the (canonical) surface S 
provided with the lattice defined by the graph V. We attach a two- level system (which is called a qubit in the language 
of quantum information and represents a spin-1/2 particle) to each edge e of the graph. Label its basis states as |± e ). 
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Consider the following Hamiltonian which is a sum of local operators attached to the vertices and the faces of the 
lattice: 

H = -E®^ e) -E®^ e) > ( n ) 

v e9u / eg/ 

with a x |±) = |=p) and a z |±) = ±|±). Calling the operators A„ = ® e ^ v Ox and B/ = <X> ee /eri e ', we first notice that 
all these operators commute with each other. We can thus diagonalize them simultaneously. Ground states \ipo) are 
then states which diagonalize all A t , and B f with the highest eigenvalue: 

A„^o) = B/I^o) = \i> ). (12) 

As shown in [5[ (and reviewed in [(|), this spin system is equivalent to BF theory based on the discrete gauge group 
7,2- Generalizing this to systems with a higher number of levels attached to each edge allows to reformulate in a 
similar fashion BF theory for an arbitrary discrete group. 

Ground states correspond to physical states in BF theory. The holonomy can either be + or — along an edge. 
Then the A v \tJjq) — IV'o) condition implements the Gauss law imposing gauge invariancc at each vertex, while the 
face condition B/|?/>o) = |"0o) imposes the flatness of the Z2 holonomy. Thus the Hilbert space of ground states for 
an orientable surface of genus g has a 2 2g degeneracy, with the holonomy around all 2g cycles of the surfaces being 
free. 

Finally, using the stabilizer space methods, the entropy of one (bounded connected) region of the lattice was 
computed in any ground state [6[. It was shown that it docs not depend on the particular ground state (i.e non- 
contractible cycles do not matter) and it is simply related to the perimeter of the region's boundary as 

S = n L - 1, 

where is the number of spins in the perimeter of a region (for a 2d square lattice). 

In the present work, we generalize these entanglement calculation to the case of the continuous Lie group, making 
explicit calculations in the case of the group G = SU(2). We do not use the same methods as developed in the analysis 
of the spin systems but exploit simple loop quantum gravity tools. 



II. ENTANGLEMENT: GENERIC STRUCTURES AND ENTROPY-BOUNDARY LAW 
A. The setting: flat spin networks and choice of independent loops 

Let us consider a fixed connected oriented (abstract) graph T and a spin network wave functional ip(g e ) living on it. 
We would like to study the completely flat wave functional, imposing that the holonomy along all loops of the graph. 
This corresponds to the flat physical state for BF theory for a trivial topology (a two-sphere) of the canonical surface 
S. It also gives the physical state for a non-trivial topology but with the additional constraint of trivial holonomies 
around every cycle of the canonical surface. To truly consider non-trivial topologies, we would need to alow for 
non-trivial holonomies around some of the loops of the graph (which would then represent the non-contractiblc cycles 
of the canonical surface). 

In order to study the entanglement properties of this completely flat state, we write the flat wave functional ipo(g e ) 
as a product of (5-functions on the group SU(2) for every loop of the graph. However, such a naive product would 
obviously lead to infinities due to redundant (5-functions and we need to consider the product over independent loops. 
We formalize this using by the gauge fixing procedure for spin networks that was developed in [ill ]. 

Considering a generic spin network functional VKffe) on T, we can gauge fix the SU(2) gauge transformations ((H) 
acting at the vertices by introducing a maximal tree T on the graph. T is a connected set of edges on T touching all 
vertices of T but never forming any loop. If we call E and V respectively the number of edges and vertices of T, then 
the number of edges in T is exactly V — 1. On the other hand, the number of edges not in T is exactly the number of 
independent loops on the graph L = E — V + 1. The gauge fixing consists in fixing all the group elements g e living 
on edges e £ T belonging to the tree to the identity. More precisely, we choose a reference vertex vo. Then for all 
vertices v, there exists a unique path [vq — *■ v] linking vq to v along the tree T. We perform a gauge transformation 
with parameters: 

h v = g e . 
ee[«o— ►«] 
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For all edges on the tree e <G T, the resulting group elements h s ^g e h t ^ are set to the identity 1. For all edges not 
belonging to the tree e ^ T, this defines a loop variable G e = h s ^g e hZ e y which is the oriented product of the loop 
C e = [vq — > s(e)] U e U [<(e) — > vq] starting at Vq and going to s(e) along the tree T, and then coming back to vq along 



the tree from the the vertex t(e). The procedure ensures that (see for more details) the wave functions evaluated 
on the original <7 e 's is equal to its evaluation on the G e ^T while setting the other group elements to 1. 
Finally, we define the completely flat spin network state as: 

M9e) = I] (13) 

It is straightforward to check that the resulting state actually does not depend on the choice of neither the reference 
vertex vq nor the maximal tree T. Moreover it truly imposes the condition that the holonomy around any loop of the 
graph r is constrained to be equal to the identity 1. 

Such a distributional state is obviously not normalisable for the kinematic inner product defined with the Haar 
measure J dG e . It is indeed L 1 but not L 2 . To deal with it, we need to regularize it. The method we will use is 
the standard one when working in loop quantum gravity and spin foam models Q: expanding the state in SU(2) 
representations, we will introduce by hand a cut-off J in the representations and then study the behavior of the 
various quantities in the large spin limit J — > +oo. 

Our purpose is to consider a bounded region A of the graph T and compute the entanglement between A and the 
rest of the graph on the completely flat spin network state. This will give the entropy of A. More precisely, we choose 
a connected region A of T. We define it as a set of Vint vertices and the E- m t edges that link them (i.e any edge 
whose both source and target vertices are in A also belongs to A). We further distinguish the boundary edges, 
that connect one vertex inside A to one vertex outside, and the E ext exterior edges who do not touch the considered 
region A. We call the interior graph r; nt the graph formed by the vertices and (interior) edges of A. The exterior 
graph r oxt is defined as its complement: it consists in both exterior and boundary edges. We define the Vb boundary 
vertices which belong to both interior and exterior graphs i.e vertices in A that touch some boundary edges. Then 
the exterior graph has V ex t = V — Vi n t + Vb vertices. 

To define the entanglement between the interior and the exterior regions, we consider the reduced density matrix 
on A obtained by tracing out all the holonomics outside A from the full density matrix: 

p(ge,g e ) = tp(ge)^(ge), pint(g e &A,g e eA) = J[dge(A\iigeeA,getA)^(geeA,g e ^A)- (14) 

We are then interested in the standard measure of entanglement defined as E = — t?A Pint log pint- We would like to 
compute it on the completely flat spin network state, which a physical state for BF theory. For this purpose, we need 
to adapt the choice of the tree T used in defining the flat state ipo to the choice of the studied region A: we would 
like the definition of ipo to respect the interior/boundary/exterior structure. 

We now choose two reference vertices "Uo in the interior (there is no problem if vq itself is a boundary vertex) and 
Wq in the exterior (for the sake of simplicity, we choose wq so that it is not a boundary vertex) and two maximal trees 
3]nt j Text respectively on the interior and outer graphs. We would like to form a maximal tree T on the whole graph 
by merging the two trees. The only issue is that the interior and outer graphs, and thus the two trees, both share 
the Vb boundary vertices. Considering a boundary vertex Vb, there exists a unique path along the tree Ti nt from vq 
to «b and there also exists a unique path in the exterior along T ext from wq to v\>. If we consider the straightforward 
gluing of the two trees T; nt U T oxt , then we obtain loops as soon as there are at least two boundary vertices of the 
type [vq — ► v h — ► wq — > — ► vo\. In order to get a tree, we simply need to break these loops. For this purpose, we 
single out an arbitrary boundary vertex and we number the other boundary vertices with i = 1, .., (Vb — 1). 
We will remove one edge from each loop [vq — ► — ► wq — > v h ^ — > v ] along T int U T oxt . More precisely, for each v h l \ 
we consider the unique boundary edge in T cxt touching v£\ Finally, we define T = (T mt U T cxt ) \ {e^} and it is 
straightforward to check that T is a maximal tree on the whole graph T. Indeed there are no loops in T. Moreover 
there exists a path in T from vq to any vertex « £ T: if v is in the region A the path is simply the path [vq — > v] 

within Ti n t while if v is outside A we consider the sequence of edges [vq — ► v h °^ — * w] U [w — > v] with the first halves 
of the path in Ti nt and the second in T ext . 

We follow the previous gauge fixing procedure and define the flat state ipo according to this chosen maximal T. We 
have one holonomy loop variable per edge e ^ T not belonging to the tree T. There are three types of such edges. 
First, we identify the edges e in the region A but not belonging to the interior tree Tj nt . Second, we identify the edges 
e outside A but not belonging to the outer tree T ext . Finally, there are the boundary edges i = 1, .., (Vb — 1). For 
all these edges e ^ T, we associate the corresponding holonomies around the interior, exterior and boundary loops 
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respectively which consist in edges / going from vq to s(e) along T then along the edge e then coming back from t(e) 
to Vo along T. We call the loops C e = [vq — ► s(e)] U e U [t(e) — > i>o] and define the corresponding holonomies: 

Veer int \r int , G e =]Jg f , Veer cxt \T cxt , H e =Y[g f , Mi = l, .., (V b - l), B t = n -9/- ( 15 ) 

The flat state is then the product of (5-functions over all these loops: 

Mg e )^ n n n <w- ^ 

eer int \T int eer oxt \T cxt i=i..(vt-i) 

We insist that the whole procedure with the choice of a tree T is simply to choose a set of independent loops in order 
to impose the flatness conditions with no redundant ^-distribution. In other words, the state ?po as defined above still 
imposes that the holonomy around any loop of the whole graph T is trivial. 

Finally, a last detail is that we can cut the loops C e for exterior edges e £ r cxt \ T ext . We can have them start at 
the exterior vertex wq instead of the reference vertex vq: we define C e = [wo — > s(e)] U e U [t(e) — > wq\. This leads 
to the holonomies H~ 1 H e H where H is the oriented product of the group elements from v$ to wo along the tree T 
(going through v^). Since the (5-functions are central, replacing the H^s by H^ 1 H e H docs not change anything to 

the definition of the flat state "00 • On the other hand, the loops C e have the advantage that there only involve edges 
in the outer graph r ext (i.e not belonging to the region A). We will therefore use this prescription for the following 
entanglement calculations. 




FIG. 1: The interior and the exterior flowers are linked by Vh boundary edges. Interior edges are shown as solid lines. 

We would like to underline that the properties (entanglement,..) of the flat state tpo do not depend on the specific 
choice of tree T that we made. Indeed the choice of tree is simply a choice to gauge fixing for the spin network 
functional. The particular tree that we defined using trees in the interior and exterior regions is a convenient choice 
allowing a clear description of the "boundary loops" which are the central objects for the entanglement calculation. 



B. Computing the entanglement for the flat state 

We now compute the reduced density matrix p lrA obtained by tracing out all exterior holonomies from the flat state 
?/V The von Neumann entropy of this reduced density matrix defines the entanglement between the region A and the 
rest of the spin network. Our main result is that the entanglement scales with the size of the boundary: 

E A = S[p int ] EE - tr A p in t log Pint = (14 - 1) «( J), (17) 

where (Vt>~ 1) is the number of boundary loops between the region A and the exterior and u( J) is a unit of entanglement 
which only depends on the regulator J (representation cut-off) and does not depend on either the graph T or the 
choice of the region A. In particular, if there is a single boundary vertex, Vb = 1, the region A is totally disentangled 
from the rest of the spin network. 

Let us compute the reduced density matrix as defined by (I14p : 

Pintle, 9e) = I ' [dg efA ^ A ] H 5^9^) ]J S(G e )S(B l )S(H e ) Y[8(G e )6(Bi)5(H e ). 

e<£A e<£T e<£T 
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First, the interior loops are unaffected by the integration over exterior edges since they do not involve any edges 
outside the region A. Second, the exterior loops involve only exterior edges, thus we have the identification H e = H e . 
This produces S(l) infinities which we re-absorb in the normalisation of the reduced density matrix. Finally, the moot 
point is what happens to the boundary loops. More precisely, we have defined: 

B i = C D^ 1 D i C7\ with Ct= JJ g f , D. t = J] g f . (18) 

felvo-Hi™ ] cr int / s [w ->v^ ] cr oxt 

With this decomposition the identification SQeg^ 1 ) for all edges e € r oxt implies that the factors 6(Bi)5(Bi) can 
be re-written as 6(C^ CoCq Ci) where the holonomy C~ 1 Cq is the (oriented) product of the group elements from 

the boundary vertex to u£ along the interior tree T; nt . Therefore, the reduced density matrix reads up to a 
normalisation (later fixed by the requirement that tr^ p mt = 1): 

(VW) 

pUg ee A,g ee A) = I] 5{G e )5{G e ) J] Sp^Cad^d). (19) 

e£T int i=l 

The next step is to compute the von Neumann entropy of this density matrix, Sa = — tr^ p- m t log p m t . We have two 
types of terms: some S(G e )S(G e ) for edges not in the tree e ^ T- lnt and some 5{g{g~ 1 ) which only involve edges in the 
tree Tj nt . Since these different terms do not involve the same group elements, they can be treated separately (more 
precisely, we can do a change of variable in the integrations replacing the g e 's by the G e 's for the edges not in the 
tree T lnii and the Jacobian of the transformation is trivial due to the left and right invariance of the Haar measure) . 
The terms 6(G e )6(G e ) corresponding to the interior loops arc density matrices corresponding to pure states and thus 
have a zero entropy. The only non- vanishing contribution therefore comes from the boundary loops and gives: 

S A = (V b - l)S[S(gg-% (20) 

where S[S(gg~ 1 )] is the entropy of the to-be-normalized density matrix a(g,g) = S(gg^ 1 ). It turns out that this 
density matrix a is actually the identity matrix and its entropy is simply the log of the dimension of the Hilbcrt 
space of L 2 functions on SU(2). However, this Hilbert space has an infinite dimension and the result requires a 
regularization. 

The regularization consists in introducing a cut-off J in the representations of SU(2). We define a regularized Sj 
function in term of the SU(2) characters 1 : 

j<J 

SAgg- 1 )^ E (ai + iJxiGMT 1 ). (21) 

jeN/2 

Corresponding we consider the Hilbert space of L 2 functions on SU(2) which decompose onto matrix elements of 
the group elements involving only representations j < J. This Hilbert space Tij is spanned by the (renormalized) 



1 We could also use the q-deformation of SU(2) which provides a natural way to restrict to a finite number of representations [TEll . For 
instance, given the deformation parameter q = exp(in/J), the highest representation is j = J — 1/2. Their g-dimensions are equal to 

^ BinW + lWJ 
J smir/J 

which is the same as the evaluation of the usual SU(2) character on an angle 8 = ir/J. When J — » oo (and q — » 1), we recover 
— > (2j + 1). Hence the q-deformed version of the reduced density matrix cr(g,g) is 

3 

Its entropy (after properly normalization of the density matrix) is 

S(p q ) = logX>- 9) ] 2 = logJV(J, tt/J) ~ 31og J+ 

3 

where the factor N(J,w/J) is introduced later when dealing with topological defects. 
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(Wigner) matrix elements \/2j + 1 D J mn (g), where m and n run by integer step from —j to +j. Its dimension is thus: 

j<J 

Aj = J2 & + = 3^ + 2J )(! + J )( 3 + 4J ) ~ 8j3 /3- (22) 

ieN/2 

It is easy to check that the non-normalized density matrix Sj(gg~ 1 ) is the identity on TCj by explicitly computing its 
matrix elements 2 and therefore its entropy is S(Sj) = log Aj. Indeed, if we introduce the normalized kcts through 

(g\j,m,n) = ^J2j + lD^ig)™, (j,m,n\j' ,m' ,n') = 6jfS mm >6 nn r, 
then the normalized density matrix 

1 j ~ J 

p J = — ^2 ^2\j,m,n)(j,m,n\, (23) 

jeN/2 rn.n 

is obviously maximally mixed, hence is maximally entangled with the rest of the system [l6j . This finally proves the 
results which we announced above (a more detailed proof is presented in Appendix A): 

Result 1. The entanglement between the region A and the rest of the spin network state for the completely flat state 
-00 is: 

S A = E A [i> Q ] = (Vb-1) logAj, with logAj ~ 31ogJ + log§. (24) 

J— KXS O 

At the end of the day, it is the number of boundary vertices that is relevant for the entanglement. It could seem 
surprising since the SU(2) group elements live on the edges of the graph and therefore the degrees of freedom of the 
theory apparently live on the edges. We would then expect the entanglement to scale with the number boundary 
edges. However, the requirement of gauge invariancc implies that the physical degrees of freedom actually live on the 
(independent) loops of the graph and not simply on its edges. Then we proved that the number of (independent) 
boundary loops is directly related to the number of boundary vertices (minus one). A last comment is that we 
would still have found the same leading order behavior for the entanglement in case we had first regularize and then 
calculate the reduced state instead of the contrary as we did above 3 . 



C. Topological defects and renormalized entanglement 



We can further generalize the procedure described above to take into account topological defects. In the context 
of 3d gravity, topological defects represent point particles. For instance, a (spinlcss) particle of given mass leads to 
a non-trivial holonomy around it: we impose 5g(g) instead of 5(g) on loops around it with the angle 8 related to the 
mass (see @, H| for more details) . 

Moreover, we expect that the leading order of the entanglement does not change when including topological defects. 
Nevertheless, for any state ip we introduce the renormalized (or relative) entanglement as a regularized entropy [l4| 
of the reduced density operator of A, Ea[$\ = <5(pm t ) 

%> int ]= km S(p J mt )-S{pi int ) (25) 

J — >oc 



2 We compute, using the orthogonality of SU(2) matrix elements: 

(v / 2fcTTD^|«5 J | y/TfpLBQ = J dgdgy/(2k + l)(2j + l^iCff^G^jfoflT 1 ) = 5 jk S ac 6 h 

3 It can be understood by noting that 



fdg5(g)8(g) = 5(1) = df A, 

jSN/2 

;egration produces the same result, 

r J 

22 djd k dgxj(g)Xk(g) = djd k Sj k =A d 



jSN/2 

and rcgularization of 5-function prior to the integration produces the same result 
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which is the difference of the regularized measure of entanglement of the region A with the rest computed in the spin 
network state ip and the entanglement computed for the completely flat state ipo- When tp only includes topological 
defects, we actually expect the renormalized entanglement to converge to a finite value as the regulator J is taken to 
infinity. 

Let us first consider inserting a topological defect in an interior loop or an exterior loop 4 , i.c replacing S(G e ) or 
8{H e ) by 8g(G e ) or 5g(H e ). In the case of an exterior loop, we will get in the reduced density matrix after integration 
over the group elements living on the exterior edges a term J dH5g(H) 2 instead of J dHS(H) 2 . Such a term only 
enters the normalization of the reduced density matrix and thus does not affect the entanglement. In the case of 
an interior loop, the modified loop constraint gives a term 8g(G e )5g(G e ) in the reduced density matrix instead of 
the original S(G e )5(G e ). This still represents a pure state, thus has zero entropy. Once again, it does not affect the 
entanglement. This proves the following result: 

Result 2. Starting off from the completely flat state ipo and then adding topological defects along some loops of the 
spin network state, if no topological defect is inserted on boundary loops then the entanglement between the region A 
and the rest of the graph does not change: 

E A [ip] = 0, or equivalently E A [ip] = (V h - l)logA, 7 . (26) 

On the other hand, if we insert a topological defect on the boundary between A and the outside, or more precisely 
replace the constraint S(Bi) by Sg(Bi) on one given boundary loop, this will modify our previous entanglement calcu- 
lation. For the completely flat state the term in p mt corresponding to that loop was J dgdg 5(gg~ 1 )5(g,ig)S(gig) = 
5{gi'g~ 1 ) which gives the totally mixed state (identity density matrix on the Hubert space of L 2 functions over SU(2)). 
This now replaced by: 

f X ($) 2 

dgdgS{gg- 1 )Sg(g l g)Sg(g i g) = J dg 5g{g i g)5g(g l g~ 1 ) = ^ (2j + ! ) *J (ffigT 1 ) > ( 27 ) 

where we used the decomposition of the distribution Sg into SU(2) representations as given by eqn. (|10p . Cutting off 
the sum over representations to J, we compute the matrix elements of this reduced density matrix: 

(V2kTTD k cd \af\VW+lDi b ) = j^S 3k [^^j 2 5 ac 8 bd , (28) 

where the normalisation factor iV(J, 6) ensures that the reduced density matrix a j has a unit trace: 

iV(^) = Sx J W 2 = (3 + 4 - ,)3in 4 ^ 3 7 (3 + 4J)9 . (29) 

One can check that doing a Taylor expansion around 8 ~ 0, we recover N(J, 6 — ► 0) = Aj. Finally, computing the von 
Neuman entropy of this density matrix gives the renormalized entanglement between the region A and the outside in 
the state ipg with one topological insertion along a boundary loop: 

Result 3. Inserting a single topological defect 5g along a boundary loop between the region A and the outside leads 
to a non-zero renormalized entanglement: 



We prove in Appendix B that this expression has a finite limit for all values of 9 when the regulator J is sent to 
oo. This renormalized entanglement has a universal value when 8 is not a rational fraction of ir: 

E A (8) = -3 + log6 w -1.20824. (31) 



4 We can further introduce any topological defect on arbitrary products of these interior and exterior loops of the type Sg (G ei G e2 G e3 ..) 
and Sg(H ei H e2 H e3 ..). As shown in the detailed proof presented in Appendix A, this does not affect the entanglement at all, which 
turns out to depend only on the boundary loops. 



10 



We can also compute this value analytically 5 for some specific values of the angle 0: 

E A (n/2) = -2 + log(3/2) w -1.5945, E a (tt/3) = -2 + log2 w -1.30436, E a {-k/A) = -2 + log3- \ log 2 w -1.248. 

We see that the renormalized entropy (and entanglement) is discontinuous. This is actually a generic feature of 
infinite-dimensional systems 0, [l3| ■ It is also known that if a reasonable constraint is put on the accessible set of 
states, such as bounded mean energy, tr Hp < oo, then the entropy will be a continuous function on this set. In our 
case the entropy is continuous if the cut-off J is large but finite. This can be motivated by appealing to the geometric 
interpretation of the representation labels. Indeed, in the context of 2+1 (loop) gravity, the spin network states 
are the eigenstates of the length operator and its eigenvalues are exactly the representation labels (up to ordering 
ambiguities). Thus requiring the finite extent of the system in space would naturally impose an upper cut-off on the 
edge representation labels. 




-0.2 
-0.4 
-0.6 
-0.8 
-1 
-1.2 
-1.4 - 
-1.6 tg 



0.5 



1.5 
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FIG. 2: The regularized entanglement E J [9] — E J A (ij>g) — E A (ipo) for different values of the representation cut-off J — 10 and 



J = 50§. 



The negativity of the above result can be easily understood if we recall the the flat holonomy state is maximally 
entangled. Therefor, the states with non-trivial holonomies are less than maximally entangled, hence the negativity 
of the renormalized entropy. 

This result extends to the more general case of several topological defects inserted along different boundary loops, 
each boundary loop contributing independently to the (renormalized) entanglement between the region A and the 
rest of the spin network state. 

Finally, we are left with the possibility of a non-trivial (2d) topology, i.e non-contractible cycles in the canonical 
surface. It is easy to see that the entanglement calculations are not affected by non-contractible cycles as long as the 
region A has a trivial topology (isomorphic to the unit disk). Otherwise, if the region A contains some non-contractible 
cycles, the holonomies around them will couple to the holonomies around the cycles outside A. As an example, taking 
the case of a 2-torus, let us consider that one cycle with holonomy G is contained inside A while the dual cycle with 
holonomy H is outside. Now, flatness of the connection does not require 8{G)8(H) but the weaker condition that 
they commute GHG^ 1 H^ 1 . Then the reduced density matrix p- lnt for the region A will be given in term of integrals 
of the type: 

a(G, G) = J dHS(GHG~ 1 H~ 1 )S(GHG~ 1 H~ 1 )f(H), 

for some central function / (invariant under conjugation). The representation decomposition of such density matrix 
involves the {6j} symbol and are more complicated to analyze. We postpone this study to future investigation. 



5 For example, for 9 = n/2, Xj(i"/2) vanishes for all half-integer values of j and is equal to (— 1) J when the representation label j is an 
integer. We can then compute the renormalized entanglement using the following exact sum: 

j<J 

log(2j + 1) = - i logTT + log(2 J+1 r(J + §)). 
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III. SOME SIMPLE EXAMPLES OF ENTANGLEMENT 
A. The 6-graph 

The simplest spin network that allows a non-trivial entropy calculation is the 0-graph. It is too simple to be 
decomposed into two non-trivial inside/outside regions, but it allows to illustrate how to compute the entanglement 
between two sets of holonomies. 

1 

— +- 




3 




FIG. 3: The oriented 9-graph 

The totally flat state is (the choice of maximal tree is a single edge T = {e 3 }): 

*(ffi,8»,Ss) = S(g 1 gs 1 )S( 9 293 1 ), (32) 

where we chose the two independent loops |ei,e3] and [e2,e3|. Notice that the (5-functions also impose that the 
loop [ei,e2] carries a trivial holonomy, gig 2 = J- Constructing the (formal) density operator ((up to rcgularisa- 
tion/renormalisation) for the edge e 3 one obtains 

P(g3,h) = J dgidg 2 6(gig^ 1 )S(g 2 g^ 1 )S(g 1 g^ 1 )5{g 2 g^ 1 ) = S(g 3 g^ x ) 2 = 6(t)S(g 3 gs x ). (33) 

The expectation values are calculated according to the trace formula (O) = tr(pO)/ tvp. Hence we regularize the 
density operator by dropping the infinite multiplicative constant and truncating the remaining ^-functions in its 
expansion in term of SU(2) representations: 

1 j ~ J 1 

p J (g3,h) = ^ E + iJxjsOasi" 1 ) = ^E E Vi* + ^^MD^ig^ 1 ), (34) 

JseN/2 33 "3,113 

where the sum on js is over all half- integers and the D 3 mn (g) are the matrix elements of the (Wigncr) matrix repre- 
senting the group element g in the representation of spin j (with m and n running by integer step from —j to 
The normalization constant is determined by the trace condition tr p = 1 , so 

}<J 1 

N = Aj= Yj ( 2 i + I) 2 = g(l + 2J)(1 + J)(3 + 4J) ~ 8J 3 /3. (35) 

j'eN/2 

Introducing the following ket notation (as in [loj). 



(g\j, m, n) = yj2j + lD 3 mn (g), (j, m, n\f, mf, n') = %'<W'<W, (36) 
the state induced on the edge e 3 can be written as 

p 3 = A~E E b3m 3 n 3 )(j3TO3n 3 |, (37) 

J OS ra 3 ,n 3 

For a pure state \ty) the von Neumann entropy of is the measure of entanglement between the edge e 3 and the 
rest (see e.g [HI)- Hence 

S(*|3 : 1,2) := S{p() = logAj - 3 log J + log(8/3). (38) 

The result is obviously divergent when J — > oo. We can check that computing the entanglement between the edge e\ 
and the rest (edges e2,e3) gives the same result. Moreover, as a consistency check, we calculate the complementary 
density matrix p(g x , g x ; g 2 , g 2 ) , 



P(i2)(ffi,<?i;32,£b) = / dg^{gi, g 2 , gs)^(gi, 52,33) = S(Gi 2 )6(Gi 2 )6(g 1 g 1 ) (39) 
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The loop G12 = <72<?i 1 is the analog of an interior loops while <?i<? 3 1 plays the role of the boundary holonomy. Finally, 
as the reduced state p( 12 ) is the direct product of the pure interior state 5{G\ 2 )8{G\ 2 ) and of the mixed state S(gig^ ), 
its entropy is determined only by the latter and we see that S(pj ) = log Aj = S(p j ) as expected. 

We can further introduce a topological defect along one of the loops of the 0-graph. Let us consider the state 

*e = Se(g l9 3 1 )5(g 2 gs 1 ), (40) 
which satisfies the following constraints, for all j £ N/2, 

Xi(5iff 3 " 1 )*e(ffi,32,53) = Xj(9i9 2 1 )^e(gi,92,93) = Xj(0)#e(si,52,ff3), 



Expanding the ^-distribution in SU(2) characters, we compute the reduced density matrix for the edge e\\ 



^91) = E ifrjxAai-gi 1 )- (41) 

After truncation of the sum over representations and proper normalisation, this gives in the ket notation: 
(i) 1 \ - \ - Xj(g) ,. ... . ^t^_V- 2/ m (3 + 4 J) sinfl-sin(3 + 4J)fl 

= MjT)2^2^WTTr bmn){jmnl m9) = = I^Te ■ (42) 

j<J m,n v ' j<J 



For (9 = 0, we recover the previous reduced density matrix computed for the flat state. This leads to a regularized 
entanglement as we obtain in the generic case: 

25(*,|1 : 2,3) = lo S N(J,9) - ^(9) log ^L. (43) 

On the the hand, we can also compute the reduced density operator p\ } for the edge e^'- 

P { e } (92,~92) = [ dg 1 d g3 5e(gig^)6(g 2 g^)6g(g 1 g^)S(g 2 93 1 ) = ^{g^ 1 ) / dGS e (G) 2 . (44) 



Up to a normalisation, it is actually equal to the density matrix p^ computed above for the flat state, pf (32,32) oc 
3{929 2 1 ) — P( 2 )(92,'g 2 )- Similarly, we obtain that (33,33) oc 5{g^g^ 1 ). That means that imposing a non-trivial 
holonomy around the single loop [ei, e 2 ] by beigig^ 1 ) does not influence the entanglement for the edge e 3 vs. [ej., e 2 ] 
or for the edge e 2 vs. [ei, 63]. 



B. Further examples 




FIG. 4: Interior edges are drawn as broken lines, exterior edges are solid. A boundary loop is indicated by a thin line. 



A simple graph that is represented on Fig. UK has an equal number of all three types of edges, E ln t = E\, = E cx t = 2 
and allows a simple inside/outside distinction. The holonomies along the internal, boundary and external edges are 
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denoted by <?, q, and h, respectively. It has two internal and two external vertices, so the three independent loops 
are produced by one internal, one external, and one boundary loop. The loops are marked on Fig[4^L. The flat wave 
functional is given by (taking into account the orientation of the edges) 

*o = S(g 1 g 2 )S(h 1 h 2 )S(g 1 q 2 h 1 q 1 ). (45) 

Computing the reduced density matrix for the two interior edges gives: 

pf n l(9,9) = S(G)6(G)S(g 1 gT 1 ), (46) 

where G = g\g 2 is the holonomy along the interior loop. This interior density matrix decomposes into a direct 
product of a pure state on the interior loop and the mixed state SQjig^ 1 ). According to the previous calculations, 
the regularized entropy is simply 

S(p*°)=\ogAj. (47) 

It is easy to check that a different choice of loop or increase in their number only results in the appearance of additional 
8(1) factors that do not alter the regularized entropy. 

Let us look at the possible addition of a new boundary edge without changing the number of boundary vertices 
as in FigUjD. The new boundary edge actually does not create a new boundary loop but simply a new exterior loop 
(since it does not involve any interior edge). Thus it should not contribute to the entanglement. Indeed the new flat 
state is: 

$0, q, h) = 5(gig2)6(hih2h 3 )d(giq2hiqi)S(q3h 3 q 2 1 ). (48) 

As expected, it results in the same reduced density matrix for the interior, p 1 ^ = p 1 ^, thus also S(p 1 ^ t ) = S(p 1 ^ t ). 
Coming back to the original graph in Fig[^i, we impose a non-trivial holonomy around the boundary loop: 

*e = 5(g 1 g 2 )5{h 1 h 2 )5 e {giq 2 h 1 q 1 ). (49) 

The interior density operator is 

pf n l(g,g) = 8(G)5(G)8 e (g 1 g^), (50) 
and it leads to a renormalized entropy equals to E(ff) as we computed in the previous section (see resultf3]). 

IV. CONCLUSIONS 

We computed the entropy for a bounded region on a physical state of BF theory. Indeed, being solvable and lacking 
local degrees of freedom, BF theory allows for explicit calculations and a precise analysis of the relationship between 
boundaries and degrees of freedom. Looking at physical spin network states, we showed that the entanglement between 
the two regions for an arbitrary bipartite splitting of the spin network only depends on the structure of the boundary: 
the entropy simply scales with the size of the boundary. More precisely, we proved that the entanglement grows with 
the number of boundary vertices, E cx (Vt, — 1). We also showed that the introduction of topological defects does 
not affect this result at leading order and we computed exactly the finite entropy difference due to such particle-like 
defects on the boundary. 

Technically, we developed the necessary mathematical tools required to analyze the graph structure of a flat distri- 
butional spin network state and to regularize the resulting wave functional and entropy calculations. These tools are 
also relevant to entanglement calculations for the spin systems used for topological quantum computation such as the 
Kitaev model [!, || . 

We hope to apply these techniques to study the precise relation between gauge breaking and entropy on one hand, 
and to compute the entanglement for a more general class of spin network states relevant for loop quantum gravity 

APPENDIX A: DETAILED PROOF OF THE RESULTS 

In this Appendix we discuss the structure of independent loops of T and prove our Results 1-3 in a more general 
setting. The connected graphs T^t and r out (that consists of of the outer vertices, V ou t = V — Vi n t and E ext exterior 
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edges) contribute Ei nt — V{ n t + 1 and E cxt — V>ut + 1 independent interior and exterior loops, respectively. It is 
convenient to introduce a further classification of the remaining E D — 1 loops. A true boundary loop in addition to the 
boundary (and, possibly) exterior edges includes at least one interior edge. Loops that contain at least one boundary 
edge, but no interior ones, we (for the lack of a better term) name frontier loops. Exterior loops that were defined 
above, together with the frontier loops form all the independent loops of r cxt (as defined in Sec. IIA it consists of 
VcKt = Vout + V a exterior and boundary vertices, and E cxt + E a exterior and boundary edges). 

From the definition of frontier loops it follows that their number in the set of independent loops is given by 

L f = (E ext +E h - V out - Vb + 1) - (E cxt - V out + l) = E h - V h , (Al) 

where the total number of the exterior loops is 

L ex t = -Ecxt + E D — 14xt + 1 (A2) 

Finally, the number of boundary loops in the set of independent loops of T is 

Lb = Vi, - 1. (A3) 

Fig. Q] illustrates the loop construction that is used to define the new variables. First we select L- lnt independent 
interior loops. Next, L = Vb — 1 boundary loops are constructed as follows. We select Vb different paths (vo Vi Wq) 
that link the two representative vertices, which are distinguished by the boundary vertices v\,V2, ■ ■ ■ ,vy b they are 
passing through. The numbering of boundary vertices induces a natural numbering of the boundary loops that are 
formed by combining the sequential paths (voViWo) and («o u>o). Finally, the set of L independent loops is 
completed by a necessary number of the external loops. 

The i-th boundary loop carries a holonomy Bi = giQig~^ x , where gi is the (internal) holonomy acquired on the way 
from vq to Vi, and the outer holonomy Qi is acquired along the path (ujioofj+i). The holonomies Qi are the first L a 
new exterior variables. Holonomies Hi around i = 1, ... , L ox t exterior loops (calculated with respect to wq) provide 
another L cxt independent variables. Since 

Loxt + Lb = E cxt + E D — V out < Loxt + Lb, (A4) 

we complete the set of outer variables by adjoining to the list the holonomy Qy b which is acquired by going from vv b 
to wq along the edge VoVv b wo, and picking additional variables r. 

To express an arbitrary holonomy R as a product of L independent holonomies and their inverses, all of them 
should be calculated with respect to some fixed reference point, such as vq. Hence, we pick up an arbitrary edge (say, 
Vv b ), so with respect to vq the outer flower holonomies become Hi — > gv b Qv b HiQv 9v • 

Holonomies Gi, i = 1, . . . , -L n t ; around the independent interior loops are part of the set of new internal variables. 
When Vb = Vnt, the representative interior point is on the boundary, and can be taken to be vy h . Then gy h = 1, and 
the set of interior variables is completed by the holonomies gi, i = 1, Vh — 1. Indeed, they resulting set is independent 
and 

L^ + (V b - 1) = E iat . (A5) 

On the other hand, when Vnt > Vb, we add a non-trivial gv b to the set of interior variables, and possibly more 
additional variables g. 

In these new variables the flat state ipo becomes 

-^int -^b ^cxt 

M9,h,q) = l[S(G r )l[S(g s Q s g-^)l[6(H t )l[S(Ri), (A6) 

r— 1 s—1 t—1 i 

where the last product contains the delta- functions of additional holonomies. The index i can cover all possible 
remaining loops of T. Those holonomies (calculated with respect to vo) can be expressed as 

R^igvM^ II <£ P *' r {9»Qs97+iT Pt -'H?' t -' (gvM^, (A7) 

Pi-{r,s,t} 

where number and ordering of the factors, as well as their powers 7r a = 0, ±1 are determined by the decomposition of 
Ri. The variable Ri is defined by the same equation with g — > g, G — > G. 
As a result, the interior density operator of Result 1 is 

p int (g,g,G;g,~g,G) = l[S(G r )S(G r ) / dQdHdr [] % s Q s3 ; + 1 1 )<5(.9 s Q s g; + 1 1 ) [J 8 3 (H t ) [J S(R i )6(R i ). (A8) 
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The integration results in 

Lint -C-b 

p int (g,g,G;g,~g,G) = J] S(G r )S(G r ))8(l) L ^ ]J 6(g^ l9s g s ~g s+1 ) J] S 2 (l), (A9) 

r—l s—1 i 

where the conditions g s Q s g^+\ = 1 and the factors 5{g~^ 1 g s g s g s j r \) and 5(G r ), were used to transform i?j to the 
identity. Dropping the infinite constants and making the final change of variables g' s = .g^+iffs brings the interior 
density operator to the form 

Lb 

Pi ^(g,g,G-,g,~g,G) = a int (G,G)l[S(g , s g'- 1 ), (A10) 

s=l 

which completes the proof of Result 1. 

Result 2 deals with the situation where Lb independent boundary loops that carry trivial holonomies can be 
chosen. In this case the state is given by 

-^iiit Lb L cxt 

m = n ^{Gr) n s{jg a Q.g7+i) n ti^m* ( au ) 

r s=l i=l i 

for some assignment {9} to the external and internal independent loops, and any compatible assignment to the 
dependent loops. 

After the integration over Q it reduces to 



Pint 



<J^ t } (G,G)f[6(g^ 1 g s ~g s g s+1 ) f dH% t {H t ) Y[S e M6 ei (R'i), (A12) 

s=l i 

where R[ = Yipi {r t}Gr Fz r H^ F ' does not contain any of the interior variables. Hence the reduced state 

L h 

pH2 = o$ } (G, G)f(9, G)f(6, G) J] S(g^ l9s ~g s ~g s+1 ) (A13) 



s=l 



has the same form as in Eq. (|A9|) . □ 
Finally, in the setting of Result 3 we pick the first boundary loop to have a non-trivial holonomy. Then 

Lint L-v L cxt 

M9,h,q) = Y[S(G r )6e(g 1 Q 1 gz 1 )l[5(g s Q s gjl 1 )]j6(Ht)l[5 ei (R i ), (AU) 

r=l s=2 t=l i 

where 6i is either 0,or (9, depending on the loop. Integration over all variables but Q\ results in the reduced state 
being given as 

P f nt ex a int (G, G) J] S(g^ l9s ~g s g s+1 ) f dQ l 8^ + \g l Q l g^ +1 (g l Q 1 ~g^ 1 ), (A15) 

s=2 

where n is a number of "other" loops that carry holonomies of the class 9. Since 6@(g) = 8e(g)cg, where the infinite 
constant can be expressed in a regularized fashion using the formulas of Ec. LA, the reduced state finally becomes 

L b 

P Ug,g,G;grg,G)=a h AG,G)5g(g[~g[- 1 )]JS(g'J- 1 ), (A16) 

s=2 

which trivially exhibits the desired entropy. 
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APPENDIX B: COMPUTING THE RENORMALIZED ENTANGLEMENT 

Here we prove that the renormalized entanglement E(8), as introduced in the Sec. II. C, 

EM = Jjm SOU) -Sfti.) = ^lo g «- ^g XjW *lo g ^L 

converges for any value of the angle 9. Because of the symmetry 6 <-> 7r — 9, we can restrict ourselves to < 9 < ir/2. 
Since 

N A 6 ) ~ — T2> A J ~ 8j3 /3, (Bl) 
sm 9 

it follows that the limit exists if asymptotically 

S(pf/ t ) ~ 3 log J + f{9) + 0(1/ J), (B2) 

for some function f(9). 

We have to consider two different cases, depending on whether the class angle is a rational or irrational fraction 
of 7r. In the first case the class angle is 9 = p/qit, with relatively prime p,q £ N, and there are only finitely many 
different values that sin 2 n9 takes. It is periodic with period Kg = q, and 



£sin 2 (W?) = 9/2. (B3) 



To simplify the notation we further take J = kq, k G N, and consider the limit k — > oo. 
As a result, the entropy 

1 2kq+1 sin 2 T7# 

S(p?i)~logJ-21ogsin#-- sin 2 log ^4-, (B4) 

J *— * n z sin 

n— 

becomes 

- 2fe (J-l 

5(pfi) -log J+ r ^^(2sin 2 ^)log(. 7g + - sin 2 (Z0)log[sm 2 (/#)]) + +0(1/ J), (B5) 

q j=0 1=0 

where we expressed the summation index n by n = jq + I. The entropy finally reduces to 

S(pH) ~ log J+ A|(i g[(2fe)!] + 2fclo g(7 ) + /(0) +0(1/ J) = 31ogJ + />) +0(1/J), (B6) 

which establishes the claim for 9 = p/qw,p, q £ N. 

For a generic value of 9 we establish the limit by using the Eulcr-Maclaurcn integration formula for sums. The 
asymptotic behavior of the entropy is 

S{pf n {) ~ log J +j(^J^ dnsin 2 {n9) logn - dnsin 2 (n0) log sin 2 (n6) \ + /(0) + 0(1/ J). (B7) 
Using the known integral J* dx sin 2 x log sin 2 x = tt{\ — log 2), the second term becomes 

i r 2J 

- dn sin 2 (n0) log sin 2 (n0) - 1 - 2 log 2. (B8) 
^ Jo 

The first integral has a closed form that involves sine and cosine integral functions, but the relevant part is simply 

1 r 2J 

- dn sin 2 (n0) log n 2 + 2 log 2 + 2 log J + . . . . (B9) 

./o 



hence the limit exists for all 9 and for the irrational fractions of 7r 

E A (^e) = -3 + log6. 
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